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ABSTRACT 

When scale symmetry is implemented into hidden local symmetry in low-energy strong in¬ 
teractions to arrive at a scale-invariant hidden local symmetric (HLS) theory, the scalar 
/o(500) may be interpreted as pseudo-Nambu-Goldstone (pNG) boson, i.e., dilaton, of spon¬ 
taneously broken scale invariance, joining the pseudo-scalar pNG bosons vr and the matter 
fields V = (p,uj) as relevant degrees of freedom. Implementing the skyrmion-half-skyrmion 
transition predicted at large N c in QCD at a density roughly twice the nuclear matter den¬ 
sity found in the crystal simulation of dense skyrmion matter, we determine the intrinsically 
density-dependent (IDD) “bare parameters” of the scale-invariant HLS Lagrangian matched 
to QCD at a matching scale Am- The resulting effective Lagrangian, with the parameters 
scaling with the density of the system, is applied to nuclear matter and dense baryonic matter 
relevant to massive compact stars by means of the double-decimation renormalization-group 
Viowk formalism. We satisfactorily post-diet the properties of normal nuclear matter and 
more significantly predict the EoS of dense compact-star matter that quantitatively accounts 
for the presently available data coming from both the terrestrial and space laboratories. We 
interpret the resulting structure of compact-star matter as revealing how the combination of 
hidden-scale symmetry and hidden local symmetry manifests itself in compressed baryonic 
matter. 
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1 Introduction 


In a preceding note [1J, the notion that the /o(500), the lowest scalar listed in the particle data 
booklet, is a dilaton arising from the spontaneous breaking of scale invariance in QCD (2| was 
implemented into hidden local symmetry (HLS) [3] of the light-quark vector mesons V/ = 
that embodies the non-linear realization of chiral symmetry into a scale-invariant 
hidden local symmetry theory (sHLS for short) and the resulting Lagrangian was subjected 
to the vacuum change due to the density of baryonic matter. 

In this paper, we confront the resulting formalism with what’s known of normal nuclear 
matter and make predictions on properties of dense matter appropriate for massive compact 
stars. 


Since the basic premise for the effective Lagrangian that we shall employ, sHLS, is fully 
expounded in PQ|, we shall eschew details and limit ourselves here only to what are essential for 
the calculations that we make. We shall follow closely the procedures given in [T|. The only 
issue that was not given an adequate comment in |lj is the place of an infrared (IR) fixed point 
postulated in [2] in QCD with the number of flavors Nf ~ 3 as needed in nuclear phenomena. 
It is perhaps worth making a brief remark on it. It is argued in [2] that the notion that QCD 
has an IR fixed point for Nf = 3 with the resulting “scale-chiral symmetry” solves some of the 
long-standing puzzles in particle physics that involve “light” scalar excitations. For instance, 
phrased in terms of a scale-chiral counting rule generalizing the chiral counting rule of chiral 
perturbation theory, it gives a surprisingly simple explanation of the A/ = 1/2 rule, accounts 
for the mass and width of the scalar /o(500) etc. As stressed in [dj, it could also resolve long¬ 
standing conundrums in nuclear physics involving a low-mass scalar. Unfortunately, however, 
there is, so far, no convincing proof that the three-flavor QCD has an IR fixed point: Neither 
lattice nor model-independent approaches have uncovered it. This is in contrast to QCD at 
N f 8 being studied for dilatonic Higgs model for going beyond the Standard Model (for a 
recent summary, see |4] j# 5 |. T his does not imply that an IR fixed point advocated in [2j is 
ruled out. As argued in [2j^ 6 . an IR fixed point at which scale-chiral symmetry is realized in 
Nambu-Goldstone mode has not yet been probed by the lattice work. 


In this paper, we take the point of view that in dense matter, the scale-chiral symmetry of 
the sort advocated by [2] could be present as an “emergent symmetry.” This is in some sense 
similar to hidden local symmetry which plays an equally important role in our calculation. 
The notion of hidden local symmetry which gives the famous “VD(vector dominance)” and 
KSRF relation makes sense only if the vector meson V/ is light. There are two known cases 
where the “lightness” of V/ is realized. One is the presence of the vector manifestation (VM) 
fixed point at which the vector meson mass goes to zero as does the pion mass (in the chiral 
limit) [3]. The other is supersymmertric QCD in certain parameter space [5|. In what follows, 
the VM fixed point, emerging at high density, will play a key role. In a similar vein, a scalar 
of ~ 600 MeV fluctuating around an IR fixed point will figure crucially in the equation of 
state (EoS) for dense matter. 

It is rather intriguing that the two symmetries we are combining, i.e., scalar symmetry 
and (vector) local symmetry , are hidden in baryonic matter (as in the beyond-the-Standard- 
Model regime m and seem to emerge at high density. 

^Private communication from Koichi Yamawaki 

^®Also in private communication from R. Crewther and L. Tunstall 
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2 Scale-invariant HLS Lagrangian 

The Lagrangian we will consider, sHLS, simplified from pQ, takes the form 

CsHLs(U, X ,V») » C^\ s {j-) 2 + f -f{^-fTr{MU^ + h.c.) + --- 

+ \d^xd^X + V{x) ( 1 ) 

where x = fo e<T ^ a is the “conformal compensator” field with a the non-linear dilaton field 
and f a = (x) the vev (either matter-free or in-medium), the chiral field U consists of (L,R) 
fields as U = e f* = £ l £r with ir = ■ tt, M is the quark-mass matrix representing 

chiral symmetry breaking which also breaks scale symmetry, fo a is the medium-free-vacuum 
expectation value (0|x|0) and V{x) is the dilaton potential that encodes the spontaneous 
and explicit breaking of scale invariance. For simplicity, the HLS Lagrangian is given to 
0(p 2 ), with the ellipsis standing for higher scale-chiral order terms. Note that we are taking 
the approximation q ~ 1 in the notation of Ref. [2j. The potential V(x) contains several 
unknown constants in [2] of which we do not need their specific forms for our analysis. 

2.1 Intrinsic density dependence (IDD) of “bare” parameters of sHLS 

In order to confront the Lagrangian (HD with nuclear matter and high density matter, there are 
three indispensable ingredients to consider. First, the baryon degrees of freedom have to be 
incorporated. Second, the “bare” parameters of the effective Lagrangian need to be matched 
to QCD. Third, strong correlations between nucleons, including possible phase changes, have 
to be included as one goes up in density. 

All three could in principle be handled - at least in some approximations such as large 
N c - using skyrmion description of baryons and baryonic matter [I]. Some progress has been 
made in this direction [6] but the mathematics required is still too daunting to arrive a reliable 
result. We shall therefore put baryons explicitly “by hand” in scale-chiral symmetric way. 
Let us call the baryon-field-implemented Lagrangian frsHLS for short. Since high density, 
n ~ (5 — 7)uq (where no is the nuclear matter density), going toward chiral transition is 
involved, the “bare” parameters need to have contact with QCD parameters. This will be 
done by matching the correlators of the frsHLS Lagrangian to those of QCD at an appropriate 
matching point A m lying below the chiral scale A x ~ 1 GeV, say, at about the p mass. It 
is important to note that the matching endows the “bare” parameters of EFT Lagrangian 
with dependence on the quark condensate ( qq ), the gluon condensate (G 2 ) etc. Since those 
condensates depend on the “vacuum,” they will of course depend on density which modifies 
the vacuum if the EFT Lagrangian is embedded into a medium. The crucially important 
point in our development is that the density dependence involved here is intrinsic of QCD, 
to be distinguished from the density dependence coming from (mundane) nuclear many-body 
correlations. This density dependence - that will play a key role in what follows - will be 
referred, as in p], to as “intrinsic density dependence” (IDD for short). 

2.2 Double-decimation RG procedure 

Now given the EFT Lagrangian endowed with the IDDs, nuclear dynamics is treated by 
renormalization-group decimation from the matching scale A m down to the appropriate low- 
energy scale where the processes we are interested in take place. For this, a highly versatile 
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tool is the Vi ow k strategy M. We will exploit it in this paper. A convenient - and successful 
- procedure in nuclear physics is the “double decimation” RG flow described in [9]. In fact, 
this procedure was used for the first time for the EoS for massive stars in m- In this 
paper, we will improve on it both in concept and in numerics, assuring consistency with the 
scale-chiral symmetry adopted in [T] . 

In the Vi^k framework, the double decimation consists of the first step from A m to the 
scale at which Vi ow k is obtained. The second step is to decimate to the Fermi sea around which 
fluctuations are computed to take into account multi-body correlations. This is equivalent to 
fluctuating around the Landau Fermi-liquid fixed point m- It is in doing these decimation 
calculations using Vi ow k that information from a topology change encoded in the skyrmion 
crystal treatment of dense matter enters. This transition involves no local order-parameter 
field and hence may not belong to the Ginzburg-Landau-Wilson paradigm but as will be seen, 
has a drastic impact on the EoS in compact-star matter. While the topology change that 
takes place in the skyrmion crystal is, strictly speaking, valid only in the large N c limit of 
QCD, it seems quite universal, visible already in the structure of the alpha particle with four 
nucleons da [i3]. It is thus highly plausible that such a half-skyrmion topological structure 
could be present in dense matter, say, above nuclear matter density. What is done in this 
paper is that this feature of changeover from skyrmions to half-skyrmions in the soliton 
description is translated into the bare parameters of the effective Lagrangian, in terms of 
changes in IDDs. It effectively demarcates the EFT Lagrangian into two density regimes, one 
for (I) n ^ ni /2 and the other for (II) n ^ ni/ 2 - The former (I) entails the “bare” parameters 
of the Lagrangian that carry the density dependence referred to as “IDD p atg” where the 
(pseudo-)NG bosons figure and the latter (II) “IDD matter ” in which the matter fields V u 
intervene. An interesting observation made in Higgs physics [Tj where also both hidden scale 
symmetry and hidden local symmetry enter is that the properties of (techni)vector mesons 
are scale-invariant. Intriguingly, it turns out also in dense matter that the p meson properties 
are scale-invariant, controlled by the VM fixed-point structure. 


3 Symmetry Energy, Tensor Forces and Topology 

One of the most interesting observables in dense matter is the “symmetry energy factor” S 
defined in the energy E per particle of nucleus consisting of P protons and N neutrons, i.e., 
A = P + N, 

E(n,x) = E 0 (n,0) + S(n)x 2 H-, (2) 

with x = (N — P)/A. Here n stands for baryon number density and the ellipsis stands for 
higher-power terms in x. The quadratic approximation is known to be reliable, so we focus 
on S. As is well-known, the symmetry energy is the quantity, representing neutron excess 
of the system, that plays a key role in the EoS of compact stars. It is this quantity that 
is strongly influenced by topology in the skyrmion picture, manifested through the nuclear 
tensor forces. 

As mentioned above, a robust feature of skyrmion description of nucleonic matter - that 
we shall exploit in what follows - is that there is a “changeover” from a state of skyrmions to 
a state of half-skrymions at some matter density denoted ni/ 2 ■ Its presence in the skyrmion 
framework is remarkably independent of the degrees of freedom involved and is quite insen¬ 
sitive to the parameters of the Lagrangian. It is present in the Skyrrne model with pions 
only as well as in sHLS models with the vectors and/or the dilaton [12] , Precisely at which 
density the changeover takes place is, however, model-dependent and cannot be pinned down 
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precisely in the present state of formulation. However the density at which the half-skyrmion 
appears, n\/ 2 , is found to be insensitive to the dilaton mass, the most uncertain quantity 
in the calculation. This feature is seen in the model in which skyrmions are put on crystal 
lattice IS HU. In what follows, our basic premise will be that in terms of the skyrmion picture 
justified at high density and for large N c , half-skyrmions could appear at some density in the 
vicinity of ~ 2uq . What ensues is a striking consequence on the symmetry energy 

Since the phenomenon considered is quite generic, more or less independent of the degrees 
of freedom involved, we can address the matter using the simplest model, i.e., the Skyrme 
model m that consists of two terms, the current algebra term and the Skyrme quartic term 
implemented with the conformal compensator field. It corresponds to dropping the vector 
meson fields and putting the Skyrme quartic term - which is of scale dimension 4 and hence 
scale-invariant - in place of the ellipsis in Eq. p. We expect the result to be qualitatively 
the same with the more realistic Lagrangian p. 

With skyrmions put on crystal, the easiest way to compute the symmetry energy is to 
rotationally quantize A-neutron skyrmion matter, which corresponds to calculating S from 
([2]) for x = 1 [16] . It is given by 

8 A/ ’ 

Here A/ is the isospin moment of inertia of 0(N C ) given by the space integral over the single 
cell of the hedgehog configuration Uq and the dilaton configuration. In the presence of vector 
mesons, the integral will also involve the mesons’s classical configurations. It is of the leading 
order in N c , with fluctuation corrections suppressed by 1/N C . 

The striking feature of the symmetry energy factor P turns out to be a cusp structure 
at the changeover density n ^/ 2 - which comes out at n j / 2 ~ (1.3 — 2.0)no- The numerical 
calculation of Equation ([3]) reveals that the S decreases monotonically as density increases 
toward 77 . 1/2 and then turns up and monotonically increases after ni/ 2 - 

Now it may be that the method anchored on crystal is not applicable to low-density 
matter. Furthermore nuclear matter at equilibrium density is known to be in Fermi liquid. 
Therefore one might object to applying the crystal skyrmion description not too far above the 
nuclear matter density. However it turns out that the cusp structure at a density at ~ 2no is 
not an artifact of crystal background and can be trusted. In fact what is highly nontrivial is 
that this feature can be easily reproduced by the microscopic structure of the tensor forces, 
in particular, the effect on the tensor forces of the topological change at ni/ 2 - For this, we 
use the fact that the symmetry energy is dominated by the tensor forces HZ]. We first write 
the effective Lagrangian P that implements the topology change at 771/2 i n the skyrmion 
description. To do this, we divide the density regime into two regions -R(egion) I and II - 
with the demarcation at 771 / 2 , 

R(egion) — I : 0 < 77 < 77i/ 2 , (4) 

R(egion) — II : 77 i/ 2 < 77 < n c . (5) 

As described in P we can translate the topology change into scaling (that is, IDD) prop¬ 
erties of the parameters of the Lagrangian P in the two regions. The principal parame¬ 
ters involved are the decay constants f^.a and masses of the pseudo-Nambu-Goldstone 
bosons, the coupling constants g P)U and masses m Pt0J of the hidden gauge fields etc. The 
specific parametrization that we extract from the strategy detailed in [T] and will be used 
for the Vi ow k approach presented in Section [5] is described in the next section. Here we 
make use of it in showing how the cusp in S can be understood in the given framework. In 
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this approach, one first constructs nuclear potentials in terms of the exchange of the meson 
degrees of freedom given in the Lagrangian with the scaling parameters. Apart from the 
IDDs in the Lagrangian, this is essentially what is done in nuclear chiral perturbation theory. 
Now in terms of our sHLS Lagrangian, the tensor forces Vt consist of 7 r and p exchanges, 
Vt = LjT + Vj,. The notable feature of V^’ p is that the two contributions, having the same 
radial form with different masses, come with the opposite sign. Thus the net tensor force 
involves a crucially important cancelation between the two components, which depends on 
the scaling properties of the two components. As will be seen in the next section, the details 
are a bit involved, but the qualitative feature is simple. 

As noted, the prominent feature of the net tensor force at n \/2 is the abrupt change in 
the slope. In R-I, the pion tensor is almost completely unaffected by density within the 
range of density we are considering, due to what one might interpret as the protection by 
chiral symmetry. This has been numerically confirmed up to ~ 5no- See Appendix |Aj The 
p tensor, on the other hand, gets enhanced as density increases due to the dropping of its 
mass. Since it comes with the sign opposite to the pion tensor, it cancels part of the pionic 
tensor. The net effect is then the tensor force becomes weaker as density increases. This 
tendency is in agreement with a variety of observations in nuclei, most spectacular of which 
is the long life-time for carbon-14 [T8j , i.e., the C-14 dating. (It should be mentioned that 
short-range three-body forces, present as contact interaction in chiral perturbation theory, 
could do the same suppression of the Gamow-Teller matrix element involved. As explained 
in m, however, this does NOT represent a different mechanism to that of [18]. It may be 
said that most, if not all, of the effect of the contact 3-body forces,largely responsible for the 
suppression of the Gamow-Teller matrix element, is encoded in the IDD included in [18].) The 
weakening of the net tensor force continues up to the changeover density ni/ 2 - At ni/ 2 , the 
tensor force stops decreasing, turns over and starts increasing, with the pion tensor becoming 
dominant. There are two mechanisms at work here. One is that the change of parameters 
that takes place at n 1( / 2 strongly suppresses the overall strength of the p tensor force although 
the mass continues dropping. The other is that in R-II, the candidate order parameter for 
chiral symmetry is a four-quark condensate with the bilinear quark condensate suppressed (it 
goes to zero at n j / 2 in the skyrmion crystal). And the four-quark condensate is found to be 
strongly suppressed in Region II [6], Interpreted in terms of a GMOR relation for in-medium 
pion, this would imply, since the in-medium decay constant remains more or less un-scaling 
in density in R-II, that the pion mass must then decrease. As a consequence, the pion tensor 
must become stronger, even further enhanced over and above the free-space value. This will 
facilitate pion condensation, as expected at high density in crystal form. Given the abrupt 
change in the tensor force at ni/ 2 , the cusp structure in the symmetry energy found in the 
skymion model follows in an immediate way as explained below. 

In summary, there is a change in the slope of the symmetry energy factor S at ni/ 2 , a 
semi-classical result (in the sense of large N c effect) which is a robust feature in the framework 
of sHLS theory. How it manifests in nature requires a sophisticated treatment of many-body 
theory. What follows in this paper is a detailed analysis of this feature in the RG-implemented 
Viowk approach which takes into account high-order correlations encoded in Landau-Fermi 
liquid theory. 
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4 Topological Demarcation of Density Regimes 

4.1 Intrinsic Density Dependence (IDD) 

In this section we specify the effective sHLS Lagrangian that is endowed with the density 
dependence IDD inherited from QCD at the matching scale A m - As announced, we deal with 
two density regimes - Region I and II - when the system is embedded in medium. That there 
can be two regimes demarcated at a density above no is neither indicated by a general QCD 
argument nor by model-independent effective field theory arguments. It is however predicted 
in the skyrmion description of dense matter which is strictly valid in the large N c limit and at 
high density. We take this into account by interpreting, as described above, the demarcation 
as the changes in the density dependence of the effective Lagrangian that is applicable to the 
Viowk approach. 

The Lagrangians d6j) and (|22l) applicable in R-I and R-II, respectively, are written in 
Lorentz -invariant form. One may object to their form saying that they should actually 
take 0(3) covariant form in medium since the Lorentz symmetry is spontaneously broken. 
In fact the 0(3) covariant HLS Lagrangian was written down before [20]. However in the 
scheme we are using with the correlators, the density dependence of the bare parameters of 
the Lagrangian is in “vacuum-specific condensates.” These do not intervene in spontaneous 
breaking of Lorentz symmetry [21] . The symmetry breaking that breaks the 0(4) symmetry 
comes in RG decimations a la Vi ow k with the given Lagrangians. Furthermore in the hidden- 
scale-HLS framework, as density goes above n\/ 2 , the quark condensate, while supporting 
chiral density wave, goes to zero on average and the vector mass drops rapidly toward the 
VM fixed point. When these two phenomena take place, the Lorentz symmetry breaking 
decreases surprisingly rapidly [20]. Thus the pion velocity, for instance, approaches 1 quickly. 
We should point out that the situation is totally different in the absence of the vector meson 
with the VM fixed point [22] . The relativistic mean field approach with the Lagrangian with 
density-dependent parameters, popularly used in nuclear theory circles, is justified along this 
line of reasoning at high density. 


4.1.1 Region-I 


Consider first Region-I <[4|). This is the normal nuclear matter phase, extrapolated to density 
n 1/2 which can be described in a multitude of phenomenologically reliable models. Currently 
most popular is the chiral perturbative approach, i.e., two-flavor yPT 2 - The Vi ow k approach 
can be considered as an improved version of yPT 2 , in that one universal IDD intervenes and 
improves on the phenomenology in the vicinity of nuclear matter where data are available. 

To be specific while preserving simplicity, we write the in-medium “bare” 6sHLS La¬ 
grangian in a linearized form 


1 m* z 

Cj = N[i 1 ^ + igV^)-m* N + g <J a]N--V 2 v + ^-V 2 


1 


+ 7T (M- 


2 

2 m a 2 


2 V 2, + 7 ^ 7 ? • d^TT - + £* 


+ 


( 6 ) 


where the ellipsis stands for possible terms that are of higher order in chiral-scale counting 
and of higher fields and V fl = r • p /t + cj fl assumed to be flavor-17(2) symmetric. Since the 
flavor U (2) symmetry for the vectors seems to be fairly good in the matter-free vacuum, it 
should hold also in low-density regime, i.e., R-I. (At high density in R-II, however, we will find 
that the 17(2) symmetry must break down [23].) C-nm stands for the pion-matter and pion-cr 
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couplings. The matching of the EFT Lagrangian to QCD renders the pion decay constant 
f n and the dilaton decay constant f a dependent on the QCD condensates C, i.e., ( qq ), ( G 2 ) 
etc. Since the condensates reflect the vacuum structure, in medium, the decay constants 
depend on density, which will be denoted with an asterisk, f* a . As stated, this density 
dependence is an intrinsic property of the QCD vacuum structure, to be distinguished from 
density dependence that is due to standard nuclear many-body correlations. This distinction 
arises from the strategy of matching EFT to QCD. 

Following the reasoning given in [Tj, we can relate the in-medium decay constants as 


K/hn ~ f*Jf0a = $/(n) 


(7) 


where /o-n-.a- are the decay constants in the matter-free vacuum]^ "I This follows from the 
nature of Nambu-Goldstone bosons reflecting the locking of the chiral symmetry to the scale 
symmetry, that is, IDD p tvg- The pion and dilaton decay constants depend on QCD con¬ 
densates C, the former on the quark condensate a la GMOR and the latter on both the 
quark condensate and gluon condensate mm- Since there is no lattice calculation in dense 
medium, the scaling function $ is not really known from QCD proper. For low density, one 
may resort to chiral perturbation theory XPT 2 . More pertinently - and fortunately - there 
is information from experiments where the pion decay constant is measured up to no, e.g., in 
deeply bound pionic system. One immediate consequence of IDDpArc is the d(ensity)-scaling 
of the pion mass 


mf/m^ 


$ 


1/2 


( 8 ) 


How the dilaton mass d-scales is more complicated. We will return to it later. 

As for the properties of, and coupling, to the matter fields, one needs to consider the 
IDD ma ^ er , that is, due to the matching of the vector and axial vector correlators. However 
as argued in [1] and elsewhere based on phenomenology, to the order we are considering, the 

This yields 


IDD maiter can be ignored in R-I, so we can focus only on IDD p atg effects 

m* N /mN ~ my/my = <f>/. 


#8 


(9) 


To the same approximation, the hidden gauge coupling g and the oNN coupling g a do not 
d-scale 


9*/9 ~ 91/9c 7 ~ 1 - ( 10 ) 

On the contrary, the pion-NN coupling in gn-NN^^f ■ d-scale 

9ttNn/ 9nNN ~ T/. (11) 

Given the “c ~ 1 approximation” made in locking chiral symmetry to scale symmetry [I] , we use approx¬ 
imate equality instead of equality. 

^ s The matching of the vector and axial vector correlators does make the pion decay constant f n inherit the 
quark and gluon condensates from QCD but their effects are negligible. It cannot account for vanishing pion 
decay constant as the quark condensate is dialled to zero. It requires a subtle role of quadratic divergence in 
the pion loops in RG decimation. Furthermore f 0 - that locks scale symmetry to chiral symmetry - cannot 
enter into the correlators of the isovector currents we have for IDD ma tt er . On the contrary, we will see later 
the situation is entirely different in Region-II. 

#9 It should be noted that the conformal compensator trick used in this paper works differently between 
the linear 7r-nucleon coupling which is used for (US and the nonlinear coupling that figures in the fcsHLS 
Lagrangian. In the latter, the axial coupling constant gA will not d-scale and hence neither will g^NN- This 
has to do with a well-known problem of the so-called “quenching of gA in nuclei” that comes from the role 
of the A resonance in the baryon sector. The gA obtained in the linear coupling accounts for the role of the 
A that is integrated out from the 6sHLS Lagrangian. This is an old story that dates back to 1974 with the 
quenching of <?a by a A-liole mechanism. See, e.g., [24]. In this paper we will use this scaling which is not 
properly included in the IDD p atg but is required for consistency. 


c #9 
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This implies, by the low-energy theorem known as Goldberger-Treiman relation, 

9A/g A **i- ( 12 ) 

Finally we turn to the dilaton mass m*. As discussed at length in [1], the dialton being a 
pseudo-Goldstone scalar with explicit scale-symmetry breaking due to an intricate interplay, 
un-understood yet, of the trace anomaly and the current quark mass, we are unable to 
determine with confidence the d-scaling of the dilaton mass with the dilaton potential of [2]. 
If however one took the dilaton potential of the Coleman-Weinberg-type log potential just to 
have an idea, one would obtain 


m* a /m a ~ <F/. 


(13) 


In R-I, a reasonable parametrization that we shall use is 


$i(«) ~ T—-7—■ 

1 + cj n/riQ 


(14) 


The value of cj > 0 used in numerical analysis will be given in Section [6] 

Though not highly rigorous, this is supported up to nuclear matter density [25] by IDD- 
implemented Walecka-type mean-field, so we will assume it in the numerical analysis given 
below. This completely determines the bare Lagrangian (JUD • Only one d-scaling function <h/ 
is to be determined and this can be done by resorting to pionic nuclear systems and/or chiral 
perturbation theory. For quantitatively accurate agreement with Nature, however, a small 
fine-tuning on cj will be required in Section [6] 


4.1.2 Region-II 

In this region, there is no guidance either from experimental data or from trustful theory - 
except for the hidden local symmetry prediction given below in (|15f) . This makes a precise de¬ 
termination of the effective Lagrangian problematic. Thus our approach is highly exploratory 
and uncertain. What is clear is that the density dependence of the parameters must undergo 
drastic modifications as the system goes across the changeover point n i/ 2 : First chiral pertur¬ 
bation theory, formulated to work well up to nuclear matter density, most likely breaks down 
at some high density in Region-II. This is because chiral perturbation theory makes sense in 
small-fci? expansion whereas Fermi liquid fixed point approach relies on small 1/kp expan¬ 
sion [26]. Second the local U( 2) symmetry assumed in R-I is likely to break down. Third, 
most significantly, in hidden local symmetry for the p meson which would be more justified 
as the vector-meson mass drops to the level of pNG bosons, there is the vector manifestation 
(VM) of hidden local symmetry, at the approach to which the mass d-scale to zero as 

m p ~ g p ~ (qq) ->• 0 (15) 

as (qq) —>• 0, where we define g p as the hidden local gauge coupling for p to distinguish it 
from g u for oj. What is significant in this behavior is that it is the hidden gauge coupling g p 
- which is un-scaling in R-I unaffected by IDD p wg - that plays an important role. Similarly 
the pion decay constant vanishing only very near the density at which chiral symmetry is 
restored, hence in R-II approaching the density that drives the system to the VM fixed point, 
is intricately connected to the matching process [3]. This means that IDD ma « er must become 
operative in R-II; (]T]) from a phenomenological point of view, were the parameters of Region I 
to continue to higher density much above no, then the symmetry energy factor would become 
“supersoft” at a density n (3 —4)no which would require modification to gravity theory [271 . 
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There is also a possibility that the Fermi-liquid structure, assumed to hold in R-I, breaks 
down in R-II. This possibility will not be considered in this paper. 

To account for a rapid changeover at ni/ 2 in sHLS, we take the d-scaling for the p vector 
meson in R-II to be consistent with the VM 


m* p /m p ccg*/g p = § p n . 

Approaching the VM fixed point, we take the linear density scaling 


(16) 


®n( n ) ~ (! - c n n / n o) (17) 

with Cjj will be fixed to give the chiral restoration density, for rough estimate, n c ~ (6 — 7)uq. 

The density n 1 / 2 ~ 2uq may be a bit too far from the VM fixed point for this d-scaling 
(|16li to be quantitatively accurate, but one can take this as expanding around the VM fixed 
point as was done for kaon condensation that takes place at n ~ 3no [28]. An approximately 
same critical density - near n jy 2 - is arrived at by expanding around equilibrium nuclear 
matter treated as the Fermi-liquid fixed point |25j . 

If the local 17(2) symmetry for ( p,oj ) were good in R-II as it seems to be in R-I, one 
could use the same reasoning given above for p. However there is nothing to indicate that 
the symmetry would hold there. For instance, the reasoning that goes into the VM fixed 
point for the p based on correlators as given in [3j does not apply to the u meson. In fact if 
one assumes U (2) symmetry and let the oj behave in the same way as the p in R-II with the 
VM property (|17l) . both symmetric nuclear matter and neutron matter become unstable just 
above the demarcation density ny/ 2 ■ This feature, shown in Appendix [C] is the first clear 
message within the framework developed in this paper that U (2) symmetry could be badly 
broken at high density. We shall therefore relinquish the U (2) hidden local symmetry for the 
vector mesons and treat the p in 517(2) HLS and the oj in 17(1) HLS as in [23] . 

The uj mass formula takes the same Higgsed mass as that of p , 


2 j2 2 
J UJ 9 UJ > 

where m 2 p = f p g p and f u is the 17(1) analog to f p 
we define a u as 


pjii. 


#10 


(18) 

In analogy to the case of p, 


u 


yj Qbjfir- 


(19) 


Now we do not know ho w m, ., scales apart from the IDD p tvg factor f*. In fact, neither f w nor 


g u is known in medium 


#n 


In what follows in confronting Nature, we will rely on Nature 


to guide us in arriving at the properties of oj at high density. 


For the d-scaling of other quantities, we again resort to qualitative features found in the 
skyrmion crystal simulation focusing on the skyrmion-half-skyrmion changeover [6|. They 
are 


• In-medium nucleon mass m* N goes like f* which is consistent with the large N c prop¬ 
erty m* N ~ ef* where e ~ 0(N^ 2 ) is related to the scale-invariant Skyrme term in 
the Skyrme Lagrangian, hence non-d-scaling. Somewhat surprisingly, the pion decay 
constant remains roughly non-d-scaling after n-i/ 2 until very near the chiral restoration 
point. Therefore we think it reasonable to take 

m* N /m N « m*/m a « f*/f 0(T « /*// 07r « k (20) 

^ 10 Phenomenologically, a p is determined to be ~ 2.1 in the matter-free space[3]. 

^ n Note in the vacuum the near mass degeneracy of p and w gives the hint that g p « g u and a p sa a^. 
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where k < 1 is more or less non-d-scaling constant up to near the chiral restoration 
density at which it could drop to zero. 

• If one assumes that the Goldberger-Treiman-like relation with the dilaton holds, i.e., 
m * N ~ g* a fa DP, then it is a good approximation to take 

g*/g a ~ constant ~ 1. (21) 

This completes the density dependence of the Lagrangian in R-II, 

1 f i t * 

Cu = N[i 7 ^ + ig*V> l )-m* N + g a a}N--V^ + -¥-V 2 

+ ^(d^a) 2 - + - ^m* n 2 TT 2 + C nm -\ - ( 22 ) 

If U( 2) symmetry held for the vector mesons, there would be only two parameters in Region 
II, the constant n and the d-scaling factor <1 ?jj that should go to zero as the VM fixed point 
is approached. Recall that in Region I, there is only one d-scaling function 4>/. As shown 
in Appendix [Cl the symmetry is broken in medium in R-II, hence in (1221) . In this case, one 
expects two additional scaling parameters for the u, i.e., g* and a*. In the analysis made 
below, we will see how these parameters are constrained by the EoS for compact stars. 


4.2 Effect on the Tensor Forces 


In this subsection, we describe the structure of the tensor forces affected by the IDDs given 
in the previous subsection. Here the u meson turns out to affect little the symmetry energy 
factor S. To see the qualitative feature of the tensor force in medium, we use the non- 
relativistic 1) form of the tensor potential, valid in the region we are considering as 

m N 

the in-medium nucleon mass stays heavy. The tensor potential [18) 122] is given by 


Vm ( r ) 


c * 2 


t 2 S l2 I{m* M r) 


47T 




_( m M r ) 3 (rn* M r) 


+ 


3m* M r 



where M = it, p, S p ^ = +1(—1) and 


c o Oh ' O ( ff 2 • r) 

>-’12 — o -5- 


— C?1 • <7 2 


(23) 

(24) 


(25) 


with the Pauli matrices r* and <7* for the isospin and spin of the nucleons with i = 1,2,3. 
The strength scales as 

— fNM rj 9 *mnn m *M ^ 6 ) 

fNM gMNN m* N mM 

where gMNN are the effective meson-nucleon couplings. Their relations to the coupling con¬ 
stants that figure in the Lagrangians will be specified below. What is significant in (I23|) is that 
given the same radial dependence, the two forces (through the pion and p meson exchanges) 
come with an opposite sign and this well-known fact plays a crucial role. 
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First, we discuss the d-scalings of the tensor forces in medium given by IDDs and predict 
how the net tensor force scales in density. For the ir tensor force, applying IDD p 7VG(i n R-I) 
and IDD ma tt er (dominantly in R-II) to n and N, we have from Eqs. ([9]), (1111) and (1201) 


hence 




9nNN m N m w 
9nNN m N TH-jx 

x dJj 1 ( 

.-i 


K X K 


for R-I 
for R-II 


RtI 


—— in both R-I and R-II. 

m n 


(27) 

(28) 


(29) 


Thus the 7r-tensor force principally depends only on the d-scaling of m£. It turns out as 
has been assumed since a long time that the pion tensor is insensitive to density: Within 
R-I, to the extent that the small pseudo-NG pion mass is in some sense protected by chiral 
symmetry, we expect the d-scaling of R n 2 to be small. And so will be Z(m*r). In addition 
the product of the former - decreasing - and the latter — increasing - largely cancels out. 
Thus the pion tensor does not d-scale in R-I. As for R-II, the situation is somewhat more 
involved. While R n is still expected to scale proportionally to the in-medium pion mass, the 
pion mass will not d-scale proportionally to y 7 (qq) since the bilinear quark condensate tends 
to zero for n ^ n\/ 2 - The chiral symmetry is still spontaneously broken in R-II, hence we 
expect the GMOR relation, expected to hold in the tree (or mean-field) order medium, to be 
modified to 


f* 2 m* n 2 = m q (qq) + ^ c n ((qq) n ) 

(30) 

n> 1 


K 2 fL m l 2 = c ™((<7<7) n ) • 

(31) 


n>l 


We have indicated by the multiquark (or higher dimension filed) condensates possible non¬ 
vanishing contribution to the GMOR mass formula in which the quark condensate (qq) is 
vanishing. There are several proposals for specific form for the order parameter(s) |30l . 
Whatever the precise form may be, multiquark condensates are expected to be quite sup¬ 
pressed in R-II as shown in |BJ. This implies that the pion mass must decrease rapidly in R-II. 
Despite the rapid decrease of the pion mass in R-II, the pion tensor remains non-d-scaling. 
This is shown in Fig. fT2l 

Now we turn to the d-scaling of the /3-tensor force. At the mean-field (or tree) order in 
sHLS, the /3-meson mass will satisfy the KSRF formula with IDD parameters 

rn p = yfa* p g* /*. (32) 

In the matter-free vacuum, the KSRF is a low-energy theorem proven to hold to all loop- 
orders in HLS [311 [3]. Whether it also holds in medium to all loop-orders or not has not been 
proven. It seems however reasonable to assume that for a given density, this does hold with 
g* replaced by the effective pNN coupling constant g p NN- It has been shown in [32] that 

<f P NN = f; 9 ; ( 33 ) 

with F* that goes to zero at the dilaton-limit fixed point (DLFP), possibly identical to the 
IR fixed point of [2], independently of how g* d-scales. In our application in Section [5l the 
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effect of F* could in principle be included. Therefore we shall leave it out in what follows 
in our discussion, setting F* = 1. with the possibility in mind that F* effect could further 
speed up the dropping in R p given below. 

While based on the d-scaling argument in [I] and relying on phenomenological observa¬ 
tions mm, possible IDD matter effect in g* was ignored in R-I, the lDD ma tter, as argued in 
Section [4.1.21 cannot be ignored in R-II. Noting that 

m*/m p 


9j\ (f_l 

9pJ \f* 

for R-I 

<F'jj x K for R-II 


(34) 

(35) 


we find the d-scaling of R p to be of the form 


R p 


9* p nn rriN m * P 

SpNN m* N m p 



j 1 for R-I 
\ ($ p n f for R-II 


(36) 

(37) 

(38) 


What is crucially important for the structure of the p tensor force is the factor R p which 
changes discontinuously from R-I to R-II across ni/ 2 with the topology change. In R-I, 
experimentally controlled nuclear processes indicate how d-scales. It is a slow decrease 
reaching ^ 0.8 at n ~ no- On the other hand, is totally unknown. It is given neither 
theoretically nor phenomenologically. The only constraint based on HLS [3] is the vector 
manifestation fixed point at which <& p n should approach 0 in the chiral limit. If the vector- 
manifestation (VM) fixed point is ~ (6 — 7)no - which is not too far from the density of the 
interior of ~ 2 solar-mass stars - then should drop more rapidly in R-II than in R-I. 
This point was already emphasized in [9]. One can see from (13811 that there will be a rapid 
suppression of the p tensor force at n-i/ 2 • This feature is shown in Fig. [0 Just for illustration, 
we have taken <f>/ = = 1 — 0.15^ and re = 1. In the application to the EoS for nuclear 

matter and compact-star matter, more realistic d-scaling of the parameters involved will be 
used. 

It is now easy to see how the cusp structure in the S factor at ni /2 arises. In the density 
regime in the vicinity of nuclear matter, the symmetry energy factor S, dominated by the 
tensor forces, can be reliably approximated by the closure formula [33] . 

S«H(R t 2 ) (39) 

where E ~ 200MeV is the average energy typical of the tensor force excitation and Vr is 
the radial part of the net tensor force defined in Fig. [TJ The tensor force strength decreasing 
as density approaches ni/ 2 from below and increasing after above rii/ 2 reproduces the cusp 
structure [33] . We suggest that this feature provides, albeit indirect, support to the scaling 
properties formulated in a general term in [Tj. 

In Section [6] the tensor force structure obtained above, together with the d-scaling prop¬ 
erties in Regions I and II, will be confronted with the EoS of nuclear matter and neutron-star 
matter. For this, the renormalization group implemented Vi ow k technique will be employed. 
This will be briefly reviewed in Section [5] 
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Figure 1: Vr (r) = Vp (r) (ti T 2 S'12) 1 - For illustration, we take n 1/2 = 2 rio, 3>j = <& p n = 
1 — 0.15— and k = 1. 


5 Renormalization Group with V/ ow fc 


As stated in Introduction, the scale-invariant HLS Lagrangian with baryons (6sHLS) can 
be applied to many-nucleon systems in either RMF involving single decimation or more 
microscopically with double decimations involving Vi ow k- Here we briefly review the latter 
approach that will be used in Section [6j Details are found in the review articles referred in 
this paper. We will essentially follow the strategy used in CHI- 

One (in principle) starts with the NN potential Vnn gotten from the 6sHLS Lagrangian 
([6]) and ( 1221 ) with the proper IDDs. One then arrives at Vi ow k a la renormalization group 
technique[71|8j. This consists of decimating the high momentum components from the match¬ 
ing scale or A e ff to a model-space momentum scale A which is not far from the Fermi mo¬ 
mentum kF • In terms of T-matrix, this amounts to computing Viowk as 


T(k',k,k 2 ) = V NN {k',k) + -V I 

71 Jo 


Vnn (k',q)T(q,k,k 2 ) 2 

-fcw- qiq ' 


^ lowk 


(■ k' 1 k,k 2 ) = V lowk (k',k) + -V [ 

77 Jo 


klowk^ > (?)-^lowk(i7! k,k ) 


k 2 — q 2 


q dq, 


(40) 


(41) 


T{k', k, k 2 ) = T lowk {k', k, k 2 ); (A/, k) < A. (42) 

Here V denotes principal-value integration and the intermediate state momentum q is inte¬ 
grated from 0 to 00 for the whole-space T and from 0 to A for 7i ow k. 

With the given “bare” effective Lagrangian, if one wishes, one can do a sophisticated 
effective field theory calculation (such as chiral perturbation theory) to obtain Vnn■ This 
should be feasible starting with the effective Lagrangian we are dealing with, i.e., 6sHLS. For 
the exploratory work we are doing here, however, a rigorous EFT calculation is unnecessary. 
In the present work, as in pH] , we choose the Vnn to be the realistic BonnS [35] NN interaction 
with the IDD dependence encoded in the “bare” parameters taken into account. We shall 
adopt the vacuum (free-space) parameters chosen in |35| without adjustments. This is a 
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phenomenologically powerful approach, fit to experimental data in free space as well as in 
medium to the momentum/energy scale defined by the cutoff A. Because we shall calculate 
the EoS, in particular, the nuclear symmetry energy, E sym (n ) up to n ~ 5no, we shall 
use A = 3 fm” 1 [36] • The Vi ow k so obtained preserves the low-energy phase shifts in the 
vacuum (up to energy A 2 ) and the deuteron binding energy of Vnn- (For example, the 
deuteron binding energy given by Viowk of A = 2.0 and 3.0 fm -1 are both -2.226 MeV.) By 
construction, Vi ow k is a smooth ‘tamed’ potential which is suitable for being used directly in 
many-body calculations. 

The first step in the procedure is to verify if in R-I the above Vi ow k interaction can 
satisfactorily reproduce the empirical nuclear matter saturation properties (saturation density 
no — 0.16 fm -3 and average energy per nucleon Eq/A ~ —16 MeV at saturation). To do 
this, we shall calculate no and Eq/A using a low-momentum ring-diagram approach [5711381 
[39] [40, 36], where the pphh ring diagrams are summed to all orders within a model space of 
the cutoff A. 

We now briefly describe the above ring-diagram method. The ground state energy shift is 
defined as AEq = Eq — Eq 1 ' 66 where Eq is the true ground-state energy and the corresponding 
quantity for the non-interacting system Eq cc . In the present work, we consider AEq as given 
by the all-order sum of the pphh ring diagrams as shown in (b), (c) and (d) of Fig. [2] they 
being respectively a (1st-, fourth- and eighth-order) such ring diagram. In our ring-diagram 



Figure 2: Diagrams included in the pphh ring-diagram summation for the ground state energy 
of nuclear matter. Included are self-energy insertions on the single-particle propagator as 
indicated by (a), and pphh ring diagrams by (b,c,d). 


calculations, we also include HF single-paticle insertions to all orders as illustrated by (a) 
of the figure. Note that each vertex of the diagrams is a Vi ow k interaction obtained from 
a density-scaled Vnn potential. We include in general three types of ring diagrams, the 
proton-proton, neutron-neutron and proton-neutron ones. The proton and neutron Fermi 
momenta are, respectively, kp p = (3-7T 2 n p ) 1 / 3 and kp n = (3vr 2 n n ) 1 / 3 , where n p and n n denote 
respectively the proton- and neutron-density. The asymmetric parameter is a = (n n — 
Tip)/(jin + Tip). With such ring diagrams summed to all orders, we have 

A E(n, a) = f d\ E E Y m (ij, A) x Y£(kl, \){ij\V lowk \kl), (43) 

m ijkl<A 

where the transition amplitudes Y are obtaind from a pphh RPA equation [37l [38]. Note that 
A is a strength parameter, integrated from 0 to 1. The above ring-diagram method reduces to 
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the usual HF method if only the first-order ring diagram is included. In this case, the above 
energy shift becomes A E(n,a)HF = ^ E n i n j(*ilMowk|k/) where rik=( 1,0) if k(<,>)kpp for 
proton and nfc=(l,0) if k(<. >)kFn for neutron. 

The above Vi ow k ring-diagram framework has been applied to symmetric and asymmetric 
nuclear matter [371138] and to the nuclear symmetry energy [36] . This framework has also been 
tested by applying it to dilute cold neutron matter in the limit that the 1 So scattering length 
of the underlying interaction approaches infinity [39l HO ] . This limit - which is a conformal 
fixed point - is usually referred to as the unitary limit, and the corresponding potentials 
are the unitarity potentials. For many-body systems at this limit, the ratio £ = Eo/E^' ec 
is expected to be a universal constant of value ~ 0.44. ( Eq and EQ ,ce have been defined 
earlier.) The above ring-diagram method has been used to calculate neutron matter using 
several very different unitarity potentials (a unitarity CDBonn potential obtained by tuning 
its meson parameters, and several square-well unitarity potentials) [39, |T0]. The £ ratios 
given by our calculations for all these different unitarity potentials are all close to 0.44, 
in good agreement with the Quantum-Monte-Carlo results (see [40] and references quoted 
therein). In fact our ring-diagram results for £ are significantly better than those given by 
HF and BHF (Brueckner HF) [391140] , It is desirable that the above unitary calculations have 
provided satisfactory results, supporting the reliability of our Vi ow k ring-diagram framework 
for calculating the nuclear matter EoSs. 

One should recognize that the above many-body approach is essentially equivalent to 
doing Landau Fermi-liquid fixed point theory with quasiparticle correlations on top of the 
Fermi sea treated with Vi ow k with IDDs implemented, as discussed in HU. This procedure is a 
microscopic improvement on the relativistic mean-field treatment involving single decimation 
of 6sHLS Lagrangian. In the application to denser regime going into R-II, the above procedure 
will be simply extrapolated. It is most likely a valid procedure if the Fermi-liquid structure 
holds in R-II. 


6 EoS for Nuclear Matter and Compact-Star Matter 

In this section, we shall extrapolate the treatment presented above, verified up to density 
no and taken to be valid up to ni/ 2 ~ 2no, to n > ni /2 appropriate to massive compact 
stars. To calculate the EoS for nuclear matter, we use the Bonn A potential[29] consistent 
with the “intrinsic density-dependent” ftsHLS Lagrangian at the leading order of scale-chiral 
counting. Here, we should note that we fix the pion exchange potential not to scale in density 
for both R-I and R-II as we argued on the basis of the pion being nearly massless Narnbu- 
Goldstone boson. This is a reasonable assumption for a qualitative account for the scaling 
of the parameters involved. As shown in [29], the central, spin-spin and spin-orbit nuclear 
forces from one-pion-exchange are weak or negligible compared with the nuclear forces from 
the other particles while the tensor force from one-pion-exchange is strong. But, as we find 
in the appendix [A] the pion-tensor is almost independent of a density. 

We start with the “bare” parameters that figure in both R-I and R-II. They are summa¬ 
rized in Table [U As stated, there is only one constant, c/, to be determined in R-I and two 
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Table 1: The d-scaling “bare” parameters of frsHLS Lagrangian ([6]) in R-I and (1221) in R-II. 

is unknown if U( 2) symmetry is broken down as described in Appendix C. How it is 
deduced is discussed in the text. 

constants, k and cjj in R-II. For the numerical work, we will adopl P^I 

c/w 0.13 -0.20, cjj » 0.15, kw 0.7 -0.8. (44) 

Before we confront Nature, we discuss how the formalism anchored solely - with no fine- 
tuning - on the three parameters of (1441) . fares. We focus on the region R-I where there is 
only one parameter cj in the adopted parametrization 4>/ = 1/(1 + cjn/no). In this region 
near no fairly well-established experimental data are available. 

First we recall that in R-I, the scaling is governed entirely by the d-scaling of /* = (x)*, 
i.e., the IDD p jvg- Although the explicit forms of the masses involved are not known in 
terms of the condensates, (qq) and ( G 2 ), we learn from CT theory [2] that the masses of 
nucleon and dilaton are dominated by (G 2 ) whereas the masses of the vector mesons are more 
crucially controlled by (qq) in particular in R-II. This means that as the condensate (qq) gets 
averaged to zero at ni/ 2 with the possible multi-quark condensates nonzero but suppressed, 
the vector-meson masses will be more strongly affected by density than the masses of nucleon 
and dilaton. Thus one expects that c/’s for p and uj are larger than c/’s for a and nucleon. 
Now suppose we simply ignore this feature and take one universal cj. The result is given in 
Appendix [Bj One sees unambiguously that neither the equilibrium density nor the binding 
energy can be gotten with only cj parameter. As we will see below, however, only a small 
adjustment within the given range of cj’s, with the above scale symmetry feature taken into 
account, can reproduce fairly well the nuclear-matter observables. This exercise demonstrates 
that Nature seems to be unreasonably fine-tuned. 

Given that a universal cj fails quite dramatically to reproduce Nature, we make, eschewing 
an extreme fine-tuning, minimal adjustments for different mesons, within the range given in 
(I44p . to calculate the ground-state energies of nuclear matter using the ring-diagram method 
with the density scaled Vi ow k interaction as described earlier. The small differences in c/ 
may be considered as 1/N C corrections in different channels in the “bare” parameters of the 
in-medium Lagrangian. 

^ 12 The constant ci that figures in the double decimation Vi ow k approach needs not be the same as the single¬ 
decimation value found in the calculation of the anomalous orbital gyromagnetic ratio 8gi measured in Pb mi- 
In fact, it is found to be ci « 0.28, somewhat larger than the range given in (l44l> . There is no discrepancy here. 
The latter subsumes some part of quaiparticle interactions captured in the Landau fixed-point paramters F\ 
and F 2 . 
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Table 2: The precise values for the scaling parameters in R-I and R-II. The scaling properties 
shown above are consistent with the scaling of the parameters in Table [U As stated in the 
text, the pionic parameters are taken not to scale in both R-I and R-II. 


Our results for symmetric nuclear matter and neutron matter are shown respectively in 
the left and right panels of Figs. [3l The scaling parameters employed are shown in Table [2j 

For R-I, as indicated in Fig. 03 we determine c/’s to provide a satisfactory description of 
the saturation properties of symmetric nuclear matter, giving saturation energy Eq/A=- 15.1 
MeV, the saturation density n sa t = 0.16 fm -3 and the compression modulus K = 183.2 MeV. 
The compression modulus comes out somewhat smaller than the value often quoted ^ 200 
MeV. This approach predicts a softer EoS for nuclear matter than one for neutron matter 
(given below). We will return to this matter ater. Table [2] shows that the scaling in R-I 
is consistent with the expectation that c/’s for p and w should be lager than c/’s for a and 
nucleon. It is to be noted that the inequality c^ )CT < c p,ul argued for in the context of chiral- 
scale symmetry is crucial for the fit to Nature. Thus Nature seems to exercise a fine-tuning 
that goes beyond the general framework adopted in our approach for the EoS considered for 
compact stars. The coupling constants g* and $*, associated with IDD mai t er . do not scale in 
this region. 

The scaling in R-II is our main interest. It will be shown below that the scaling in 
Table O qualitatively consistent to that of Table [T] produces the EoS for nuclear matter, 
which is compatible with Nature. For R-II, we take <b p u = 1 — 0.15(n/no) for the VM 
behavior of the p with n c « (6 — 7)no and m*/m p = g*/g p = <& p ir The oversimplified 
parametrization for m* and g* may be the cause of the most likely artificial gap in the 
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ground state energy at n = ni/ 2 - To remove this gap and get the resulting EoS within the 
empirical constraint of Danielewicz[42j. we adjust the values for Ac’s of a, u and nucleon. As 
for the u properties, we take a scaling drastically different from that of the VM behavior of 
the p, say, <4)7 = 1 — 0.053(n — ni/ 2 )/no- Given that the attraction is largely controlled by 
the dilaton exchange whose mass remains unsealing or at most slowly scaling, the repulsion 
due to the u exchange is highly constrained, so that a faster decrease of the g^/g^ cannot be 
accommodated. Likewise a substantial decrease of a* would not be allowed if one were not 

Lu 

to exceed the Danielewicz constraint. Therefore we have simply taken a* = a w . 
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Figure 4: Comparison of our calculated symmetry energies E sym (n ) with the empirical 
ones of Li et al. [43] and Tsang et al. [44]. It is worth noting that the predicted symmetry 
energy manifests a shift from soft to hard at ni / 2 reflecting the classical cusp structure in the 
skyrmion description of the topology change. 


As indicated by Eq. ([ 2 |) . the nuclear symmetry energy E sym (n) is given by Eq(u, l)/A — 
Eo(n,0)/A. (Here E sym is the same as the S factor of Eq.([2]).) In Fig. [4] we present our 
calculated symmetry energies, and compare them with their empirical values. Li et al. [43] 
have suggested an empirical relation 


E sym (n) « 31.6MeV(n/no) 7 ; 7 = 0.69— 1 . 1 , 


(45) 


for constraining the density dependence of the symmetry energy. The upper (7 = 1.1) and 
lower (7 = 0.69) constraints are also plotted in the figure, labelled as A and B respectively. 
Tsang et al. [44] proposed an empirical formula for the symmetry energy, namely 


E 


sym 


C Stk (n\ 2 ' 3 
(n) = I" [no 


n \ 7i 

2 V n 0 ) 


+ I _ 


(46) 


where C s k = 25MeV, C S:P = 35.2 MeV and 7 j ~ 0.7. This formula is also plotted in 
Fig. [4] labelled as C. As seen, our calculated E sym is slightly lower than the constraints in 
the low-density region, and is close to Li’s upper constraint in the high-density region. It 
is noteworthy that the S is relatively soft at low densities and hard at high densities, the 
changeover occurring at the crossover density This is a prediction of our theory. 

Extensive studies have been carried out by Lattimer and Lim [45] concerning the empirical 


constraints on E sym and L (defined as 3u(dE sym /du),u = n/no) at density n = no- Their re¬ 
sults deduced from a wide range of observables including nuclear masses, nuclear giant dipole 
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resonances, astrophysics and neutron skins of the Sn isotopes suggested 28 ^ E sym /MeV ^ 32 
and 40 L/MeV ^ 60. Our results for them are given in Tabled indicated by ’6sHLS’. 
Our E sym is ~ 27 MeV which is slightly smaller than the empirical value of ~ 30 MeV. Our 
L value of ~ 57 MeV is in satisfactory agreement with the empirical values of Lattimer, but 
slightly lower than Li’s lower and Tsang’s constraints. (The two entries in row 2 of the Table 
[3] are respectively the L values given by Li’s lower and upper constraints.) 


E sym /MeV 

L/MeV 


27 

57.3 

fesHLS 

31.6 

65.4-104.2 

Li 

30.1 

62.0 

Tsang 

28-32 

40-60 

Lattimer 


Table 3: Comparison of our calculated E sym and L (bsHLS) at n = no with the empirical 
values of Li et al. [43j, Tsang et al. [44] and Lattimer et al [45] . 
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Figure 5: Calculated pressure in symmetric nuclear matter (left panel) and the same in 
neutron matter (right panel) compared with the Danielewicz constraints [42] , 


It is of interest and useful to calculate the pressure-density EoS p(n) and compare it with 
the empirical constraints of Danielewicz et al. [42], This EoS is given by using 


to fit Eq/A in Fig. [3] where 


and the energy density is 


E 0 

—r = a 
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(47) 
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-e(n) 

(48) 

e(n) = 

E 0 (n) 
n[ A 

■ + m N ] 

(49) 


with mjv being the nucleon rest mass. Our result for p(n ) of symmetric nuclear matter is 
shown in the left panel of Fig. [5j Our EoS is generally in agreement with the Danielewicz 
constraint, although being rather close to the lower boundary of the constraint box. 


Our calculated p(n) for neutron matter is shown in right panel of Fig. (H and as shown 
it is generally in agreement with the Danielewicz constraints. There are two Danielewicz 
constraints in this case, one for the empirical stiff EoS (upper box) and the other for the soft 
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Figure 6: Mass-radius relation of pure neutron stars calculated from the EoS of Fig. [5] (right 
panel). Our EoS does not contain the part of EoS of low densities appropriate for the surface 
region, hence cannot account for the low-mass stars. 



Figure 7: Central densities of pure neutron stars of Fig. [6l 


one (lower box). Our p(n), being near the lower boundary of the upper (stiff) box, is mostly 
in the overlapping region allowed by both constraints. 

Our neutron-matter EoS may be considered as ‘medium stiff’. Can it support a massive 
neutron star such as the 2 Mq one recently observed by Antoniadis et al. [46]? From the 
neutron matter EoS we can calculate the properties of pure neutron stars. This is of much 
interest, and could provide a useful test of our neutron matter EoS in the high density region. 
We have done so and our results are presented below. We first calculate the pressure-energy 
EoS p(e) and then various neutron-star properties are obtained from solving the Tolman- 
Volkov-Oppenheimer (TOY) equations with the above EoS as the input. (See e.g. [38]). 


In Fig. [6] our calculated neutron-star mass-radius relation is shown. Our maximum-mass 
neutron star has mass M ~ 2.07 Mq, and radius R ~ 11.7km. With the weak equilibrium, not 
included here, taken into account, we expect that the maximum mass will come down a bit. 
This calculated mass is close to the mass 2.01 ± 0.04Mq of the recently observed relativistic 
pulsar [46]. It may be worth mentioning that the central density of our maximum-mass 
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0.8 



Figure 8: The sound velocity obtained from a fitting formula that reproduces the EoS of the 
neutron matter in Fig. [3l Near the crossover density n\j 2 = 2no, there is some discontinuity - 
most likely artificial - in velocity caused by the sharp demarcation of the parameter scaling, 
which turns out to be sensitive to the way the fitting is done. 


neutron star is merely ~ 5.6no as indicated in Fig. [71 At this density, we have found that 
our EoS is within the causual limit; for example our EoS has v soun d/c ~ 0.65 at densities 
between n/n o = 5.0 and 6.0. This is shown in Fig. [8j What is significant in this result is that 
the change in the d-scaling in the parameters of the effective Lagrangian in R-II stabilize the 
sound velocity within the causal limit. Without the topology change, the causality would be 
violated at the density reached in massive stars. 

Lattimer and Schutz have proposed an empirical relation 

. . „ 9r M km , M km , . 

I ~ (0.237 ± 0.008) Afi? 2 [1 + 4.2—— + 90( ——) 4 ] (50) 

constraining M, R and the moment of inertia I of neutron stars m- In solving the TOV 
equations, we integrate outward from the center of the neutron star till its edge where pressure 
is zero. In this process we know the matter distribution at all radii and thus can calculate 
its moment of inertia I. In Fig. (9[ we compare this I (A) with the one given by the above 
relation using our calculated M and R as inputs (B). This comparison provides a check of our 
calculated density profile of the neutron star. As seen, our results are in good agreement with 
the empirical values for neutron stars with mass less than ~ IAMq But at larger masses, 
there is significant difference between the two. As of now, we are not sure about the reason 
for this difference. In our present calculation, neutron stars are assumed to be composed 
of pure neutron matter. Maybe the above difference is related to this assumption; namely 
this assumption is possibly adequate for light neutron stars but not so for heavier ones. 
For neutron stars of mass larger than ~ IAMq, the presence of other constituents such as 
protons [38] and/or strange particles may be necessary. In the following section, we consider 
the latter possibility. 
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Figure 9: Moment of inertia of neutron stars of Fig. [6l 


7 Strangeness Problems 

We discuss in this section how the approach formulated in flavor SU (2) exploited above for 
the EoS of compact stars can be applied to strangeness in compact-star matter, namely, the 
hyperon puzzle and kaon condensation problem. In the treatment given above, the effect 
of the dilaton which has a natural habitat in flavor SU( 3) is projected into the SJJ{2) HLS 
Lagrangian. To be fully realistic in addressing strangeness, one should resort to 3-flavor 
scale-invariant HLS Lagrangian with the dilaton treated on the same footing with kaons [2jj. 
To address the EoS of compact stars with strangeness duly taken into account, the Vi ow k RG 
approach would then have to be extended to three flavors. Unfortunately such a three-flavor 
Viowk formulation is not yet available. In this section, we sketch how to address some of the 
issues involved with strangeness in the two-flavor framework developed applied in Section [6J 

7.1 Hyperon problem 

In developing the EoS in Section [6J strangeness degrees of freedom have been ignored. It is 
known that if the strangeness enters in compact-star systems, then the EoS can become too 
soft and massive stars of ~ 2 solar mass cannot be supported. For instance, this applies to 
the presence of hyperons. Simple energetic considerations suggest that hyperons should be 
present at high density in compact-star matter ^48]. The lowest-lying hyperon A, with its 
attractive interaction, is estimated to appear at matter density ~ 2no with the others possibly 
appearing at higher density. What this suggests is that the hyperons could appear at about 
the same density as the one at which the half-skyrmion phase appears in the skyrmion matter. 
If this were the case, then the prediction made in Section [6] would make no sense without the 
hyperonic degree of freedom taken into account. 

Since a Vi ow k formalism for 3-flavor is not available, we address this problem using a mean- 
field approach with the two-flavor 6sHLS Lagrangian employed in Section [6j One can think 
of this approach as a “single-decimation” RG approach as introduced in [9] in contrast to the 
double-decimation applied above. This approach was applied with success to the calculation 
of the anomalous gyromagnetic ratio in heavy nuclei 5gi m- 
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Figure 10: Tadpole diagram for self-energies for the nucleon N and the hyperon A in medium. 
The loop corresponds to the nucleon scalar density n s = (NN) for coupling to a and the 
nucleon number density n = (N^N) for coupling to uj. 


We find that with the scaling property of the “bare” parameters of the Lagrangian ob¬ 
tained above, the interactions between As and nucleons become sufficiently repulsive at a 
density n 3no so as to prevent the As from appearing in the system. Our reasoning relies 
on what we shall refer to as “Bedaque-Steiner” constraint 

In a highly sophisticated phenomenological study using a Monte Carlo simulation over 
parameters that enter in the EoS for symmetric and asymmetric nuclear matter such as the 
compression modulus K and L, symmetry energies S and S\, Bedaque and Steiner obtain 
the range of density A constrained by hydrodynamic stability of the system, that ensures 
that stars with M > 2Mq could be supported [Uj. The A is then the range of density 
beyond which the in-medium A mass becomes greater than the vacuum value. One expects 
- and it is confirmed experimentally - that the A-nucleon interaction is attractive at normal 
nuclear matter density, so As can be bound in nuclear matter. In compact star matter, as 
density increases, the chemical potential difference between neutron and proton increases and 
it can become energetically favored to have spontaneous creating of hyperons in the system. 
It turns out that this can happen when density reaches roughly twice the normal nuclear 
matter density unless the attractive interaction is cancelled by repulsive mechanisms. The 
instability generated by the presence of hyperons at that low density is the hyperon problem. 
The A then stands for the range of density at which the A-interactions must be repulsive 
enough to make the in-nredium A mass be greater than the vacuum mass. The analysis by 
Bedaque and Steiner establishes that the range of A required is 1 & A/no ^ 3. 

It is feasible, with some reasonable assumptions, to calculate the effective mass of A in 
medium using the bs HLS formalism applied above. We do this using an RMF approximation 
with the Lagrangian ([6]) and ([22]) . 

In the mean held approximation, the chemical potential for the A in medium gets contri¬ 
butions from two sources, one from the IDD p tvg hi the “bare” mass parameter m\ ~ 
and the other from the potential terms coming from A-nuclear coupling via a and uj exchanges 
as depicted in Fig. [TUI 


MA = m* A - 


9*aA9*gN 

m * 2 


n s + 


9uja9uN „ 

*2 n 

mj 1 


(51) 


where n s and n are, respectively, nucleon scalar density and nucleon number density as defined 
in Fig. M The notations for (<r, uj) coupling to A and N are self-evident. The asterisk stands 
for d-scaling parameters. 

Apart from the A coupling to the mesons, the large cancelation between the a attraction 
and the uj repulsion responsible for small binding energy for nuclear matter must take place 
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Parameters for A Mass Shift 
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Table 4: The “bare” parameter scaling for mean-field estimate of A mass shift in dense matter. 
The only scaling parameter is chosen to be c/ = 0.13 as in Section [6l The vacuum scalar 
(dilaton) mass is taken to be m a = 720 MeV so as to give ~ 600 MeV at nuclear matter 
density appropriate for RMF approach. We have taken a = g u N = 12.5 and \g a A = 
g a N = tun/ fir- The empirical values mjv = 939 MeV, m\ = 1116 MeV and = 783 MeV 
are taken from the particle data booklet. The scaling ^ ~ 1 — 0.053 ” n " /2 is taken as the 
“best fit” from the analysis in Section [6J 



0.0 0.5 1.0 1.5 2.0 2.5 


n ln Q 

Figure 11: gA—'mA vs. n/n q calculated with the scaling parameters determined in the theory 
and summarized in Tabled The demarcation density was chosen for ni /2 = (1.5, 2.0)no- The 
density at which the mass shift crosses zero corresponds to A of ^49], 


also in this case. In fact, using the standard constituent quark countin g# 13 ! , we may take 
g a A ~ | g^N and g u a ~ , the 2/3 factor accounting for the two non-strange quarks in A 

vs. 3 in nucleon. Then 
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(52) 


This indicates that the A effective mass shift gA ~ ut-a will become positive near ~ 2no as in 
the symmetric nuclear matter. 

To make a rough estimate of the A mass shift in dense medium, we take into account 
the d-scaling of the parameters in the 6sHLS Lagrangian in the mean-field calculation which 
corresponds to the “single-decimation procedure” of [Qj. In doing this, it is important to 
recognize the scaling parameter c/ in this procedure could be different, i.e., renormalized, 
from the IDD coefficient entering into the double-decimation procedure with Vi ow k employed 
in Section [6j The reason is that in the single-decimation procedure of RMF, as noted above, 

, ’ tl3 In CT theory, the dilaton is a strong mixture of the quarkonium component and the gluonium component, 
so this quark counting may not be reliable. However we do not expect it to deviate much from 1. 
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the scaling function <h/ is related to the Fermi-liquid fixed-point parameters as shown in m 
and encodes certain nonperturbative quasi-particle interactions on top of the IDD effects 
manifesting scale-chiral symmetry. We ignore this subtlety that we expect to be of higher 
order fluctuation corrections. Thus we take cj as used in the Vi ow k calculation, cj ~ 0.13. The 
scaling and the constants for the dilaton and uj used in the calculation - which are consistent 
with the property of nuclear matter treated in the mean field of the given Lagrangian - are 
summarized in Table 01 

The result is plotted in Fig. [TlJ We see that /j,\ — m a crosses zero at a density 1.5 < 
n/riQ < 2.0. The result is insensitive to the demarcation density for the regions. In fact 
what comes out in the mean field is quite easy to understand. Since m\ stops dropping 
with /* stabilizing at 2no, what matters is the interplay of the ratio (^r) 2 for the scalar and 
vector mesons - with an opposite sign - multiplied, respectively, by the scalar density n s 
and by the baryon number density n. The vector repulsion wins over the scalar attraction 
as density increases in the same way as it does in nuclear matter. Although the estimate is 
admittedly approximate - and it could be done much more realistically in the Vi ow k approach 
used above, that the BS constraint [49] is met is most likely robust. We conclude from this 
that within the formalism developed in [1J and with the prescription given in Ref. [49], the 
hyperon problem does not arise in compact stars and hence the EoS discussed in Section [6] 
with hyperon degrees of freedom ignored could stay valid. 

7.2 Kaon condensation 

In the literature, the hyperon problem is treated independently of kaon condensation. We 
believe this is incomplete if not incorrect. In fact both hyperons and kaons figure together 
in flavor SU(2>) chiral Lagrangian and should be treated on the same footing. As will be 
discussed below, to 0(N®), hyperons and condensed kaons are likely to appear at the same 
density, with the possibility that higher-order corrections in 1/N C could trigger hyperons to 
appear before condensed kaons. In the preceding section, it was suggested that hyperons may 
be ignored in the EoS. The argument developed there was quite simple. In a close parallel 
to nuclear matter at high density where the repulsion between nucleons in exchange of uj 
mesons overpowers the attraction due to scalar exchanges at densities near 2no, A-nuclear 
interactions make the effective mass of a A in medium greater than that of a A in the matter- 
free vacuum. In this section, we explore whether a related mechanism could be applied to 
avoid the “kaon condensation problem.” 

7.2.1 Callan-Klebanov skyrmion 

We first address the issue as to whether kaons condense before or after the appearance of 
hyperons. At present, to the best of our knowledge, the only way this problem can be 
addressed in a tractable approximation in consistency with the basic premise of QCD 
such as large N c - is the skyrmion description in which both kaons and hyperons can be 
treated on the same footing with a same Lagrangian. This matter was first discussed in [50] 
employing the successful Callan-Klebanov bound-state model m- In this model, anti-kaons 
K~ are bound to the SU( 2) skyrmion to yield hyperons. A highly non-trivial and surprising 
observation is that this model interpolates kaons between the chiral limit (mx 0) and the 
Isgur-Wise heavy-quark limit (mx —> oo). The model can be applied to nuclear matter by 
putting the CK skyrmions on crystal lattice. It was shown [50] that put on a crystal, the 
energy difference between the lowest-lying hyperon A and the nucleon N in medium comes 
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E* a -E* n =oj* k + 0(N~ 1 ) (53) 

where the asterisk represents medium-dependence. Note that the in-medium kaon mass is of 
N c ° in the N c counting. It is fortunate that the leading 0(N C ) term and the flavor singlet 
0(N®) Casimir energy term - which is extremely difficult to calculate - cancel out in the 
difference. 

Now, the As will appear in the compact star matter when 

He > — E^j (54) 

where He is the electron chemical potential which is equal to Hn ~ Hp i n weak equilibrium. 
On the other hand, kaons will appear by the weak process e~ —>• K~ + v e when 

He > m* K . (55) 


Therefore to the leading order in N c in QCD, hyperons and condensed kaons populate compact 
stars simultaneously. Which one appears first in the single Lagrangian description depends 
on (9(1/Ay) hyperhne corrections, namely, when the skyrmion-kaon system is rotationally 
quantized. A simple quasi-particle approximation leads to 

E* a -E* n = u;* k + ^-(c* 2 - 1) (56) 


where 11 > 0 is the moment of inertia of skyrmion rotator and c* is the in-medium hyperhne 
coefficient multiplying the effective spin operator of strangeness -1. The coefficient c is highly 
model-dependent even in the matter-free space m, so it is unknown in dense matter except 
in the large N c limit and also in the chiral limit. In either or both of these limits, c* —>• 1. 
In the matter-free space, it is found to be c 2 ~ 0.5. Although presently there is no proof, it 
seems likely that c* 2 < 1 in medium, approaching 1 from below near chiral r estor ation. If 
this is the case, that would suggest that hyperons appear before kaons condens e# 14 l and they 
ultimately join in the vicinity of chiral restoration. It is however difficult to be more precise 
on this point since the effect is at 0{N^ 1 ) and at that order there are many other corrections 
of the same order, such as higher-order nuclear correlations, that go beyond the mean-field 
order as in the Vi ow k approach of Section [6l We are therefore unable to conclude that the 
absence of hyperons a la Bedaque and Steiner precludes kaon condensation. 


7.2.2 Relativistic Mean Field with AsHLS 

If hyperons do not figure in the EoS considered in Section [6] and kaons condense only after 
hyperons appear as suggested above, does it mean that kaon condensation can also be ignored? 
In order to address this question, we need a more detailed analysis within the framework 
developed in the paper. In the absence of V low ^ ' approach, the best we can do is an RMF 
approach using the ftsHLS Lagrangian with kaons implemented as “heavy” mesons with the 
scaling parameters in the SU (2) sector fixed in Section [U A similar approach is discussed 
with dilaton treated differently from ours in [52] . 

In RMF, we can follow the argument given in [25]. In mean field, the kaon in medium 
receives mass shift by the tadpole in Fig.QUlwith the left baryon ( N , A) replaced by K~. There 
is one striking difference between the baryon and the kaon. Unlike in the case of baryon where 

^ 14 This conclusion is opposite to what was aimed at or conjectured by 1501 . 
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the scalar contribution is cancelled by the vector contribution, for the K ~, both come in with 
the same sign thanks to the G-parity and the attractions add. This immediately precludes 
the mechanism that prevented the appearance of hyperons for n ^ 2no for preventing kaon 
condensation. This would mean that kaon condensation could intervene at higher densities 
where hyperons are not present. 

What can prevent this was suggested in [53] where the authors introduce higher-order 
nuclear correlations that involve repulsions in nucleon-nucleon interactions. The mechanism 
proposed in [53] is shown to push nx to 6no, above the central density of ~ 2-solar mass 
stars. This mechanism is not captured in the mean field with our frsHLS. It could however 
be captured in a three-flavor microscopic Vi ow k approach, a project relegated to a future 
research. 

• Other considerations 

(1) Of interest is the role of the IR fixed point of scale symmetry at which the dilaton 
mass is to vanish (in the chiral limit) [2]. If the IR fixed point of scale symmetry is near n cen t, 
the possibility we have ignored in this paper, then one would have to consider the possible 
breakdown of Fermi liquid structure in the baryonic matter involved as discussed in [53] . The 
break-down of Fermi liquid structure would make this problem a whole new ball-game. Kaon 
condensation in non-Fermi liquid is a totally unknown object. 

(2) An alternative scenario is that if the kaon condensation threshold density is pushed 
up by the mechanism for the hyperon solution, then condensed kaons could be in a peace¬ 
ful coexistence with non-Fermi liquid baryons in a form similar to strong-coupling strange 
quark matter co-existing with hadronic matter with no phase transition as discussed in a 
phenomenological model [55] or 3-layer structure consisting of hadrons, condensed kaons 
and strange quarks with a kaon-condensed state playing a doorway state to strange quark 
matter [56] . 


8 Conclusion 

In this paper, we have constructed an “intrinsic density-dependent” scale-invariant hidden 
local symmetry (“6sHLS”) Lagrangian with baryons included explicitly, capturing sliding- 
vacuum properties induced by densities, and applied it to nuclear matter and compact-star 
matter. In determining the “bare” parameters of the frsHLS Lagrangian, we exploited the 
structure of baryonic matter present in the skyrmion description, namely, the skyrmion- 
half-skyrmion topology change at a density above that of the normal nuclear matter, and 
determined the IDD parameters in two regions of density R-I and R-II with the demarcation 
at the topology change density n\i 2 ~ 2no- It turns out, remarkably, that the scaling of the 
parameters of the “bare” Lagrangian with which the Vi ow k RG flow is to be performed can 
be put in the concise form ^ ^ « y v ^ ~ jy- with y v = and V = (p,u). 

Apart from the quantity yy which is controlled by IY)T> maUer , the scaling of all light-quark 
hadrons in nuclear dynamics is dictated by IDD p tvg representing the locking of chiral and 
scale symmetries. 

With no unknown parameters, the properties of nuclear matter are well described by 
the Vi ow k RG approach up to the equilibrium density no and are argued to be reliable in 
R-I up to the topology change density n X j 2 ~ 2no- In R-II, in contrast, due to the paucity 
of both theoretical and experimental input, it is found to be difficult to pin down reliably 
the parameters of the bare Lagrangian. However relying on theoretical arguments based on 
the vector manifestation property of the p vector meson and the pseudo-Nambu-Goldstone 
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nature of vr and a, we were able to fix almost all except for the y u for the w-NN coupling 
due to the apparent breakdown of U (2) symmetry at high density. Assuming that the Fermi- 
liquid structure, known to be valid in the vicinity of nuclear matter density, continues to be 
applicable in R-II up the range of densities involved in compact stars, say, ~ (5 — 6)no, we 
were able to satisfactorily confront the properties of the recently observed massive neutron 
stars. We admit that were the Fermi liquid structure broken in the density regime concerned 
- which cannot be excluded in the vicinity of the possible IR fixed point |5l], our results 
could not be trusted. 

To summarize what we have found: 

The topology change that takes place in the skyrmion description of nucleons, incorpo¬ 
rated into 6sHLS Lagrangian, has a dramatic effect in the density regime n> 2 ~ 2no on 
the EoS of compact-star matter. It affects the bare parameters of the effective Lagrangian due 
to the existence of both the VM fixed point of the p meson and the IR fixed point associated 
with the dilaton. 

• It makes the nuclear tensor forces for n ^ n 1 / 2 predominantly controlled by the pseudo- 
NG pion, with the competing p tensor strongly suppressed, and induces a shift at ni /2 
from soft, as needed in heavy-ion collisions m , to hard in the EoS, specially, the 
symmetry energy, as needed for massive neutron stars. 

• The changeover from skyrmion matter to half-skyrmion matter observed in the skyrmion- 
crystal description resembles, uncannily, the smooth transition at n ~ (2 — 3)no from 
hadronic matter to strongly-coupled quark matter recently discussed [55], [58, 59]. In 
particular the half-skyrmion phase could be identified with the quarkyonic phase of m- 

• The U (2) symmetry for (p, oj) which holds fairly well in the vacuum - and presumably in 
R-I - must be broken down at high density in R-II. Otherwise there will be inconsistency 
with the properties of the observed massive stars. 

• The topology change with the consequent IDD parameter changes makes the oj repulsion 
dominate over the a attraction in R-II. This could potentially prevent the hyperons from 
appearing at a density n 6no, thus resolving the “hyperon problem.” This could also 
be interpreted as the mechanism that accounts for the observed small binding energy 
- an order of magnitude small relative to QCD scale - for nuclei and nuclear matter, 
leading effectively to a BPS structure of baryonic matter discussed in |T], seemingly 
at odds with QCD in the large N c limit. The same mechanism could prevent kaon 
condensations in the same range of density as that of hyperons but this requires further 
studies. 

• What in our view is the most significant among our observations is the origin of proton 
mass as opposed to that of quark mass. The prediction of CT theory that the proton 
mass is dominantly gluonic, non-vanishing as the quark condensate goes to zero, and 
hence chirally invariant, is, albeit indirectly, supported by the results of this calculation. 
This suggests that the mechanism for the proton mass generation lies outside of the 
standard paradigm based on spontaneous breaking of chiral symmetry. This feature is 
supported by skyrmion crystal models as well as parity-doubling baryon models. 

Finally we should mention a fundamental issue related to the scale-chiral symmetry con¬ 
sidered and its potential generalization. In considering the scaling properties of baryons and 
mesons, we have implemented only the vector manifestation of p which brings (in the chiral 
limit) 7 r and p into a zero mass multiplet. On the other hand, the scale-chiral symmetry 
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considered in this paper implies the joining (in the chiral and scale limit) of ir and dilaton 
o' into a zero-mass multiplet. The two could correspond to the same zero-mass multiplet. 
Together with the ai, they could then constitute the multiplet figuring in Weinberg’s mended 
symmetries [60]. As discussed in [I], a possible scenario could be that ir, a, p and a\ all come 
together in a massless multiplet at the chiral-scale restoration density. In the scaling behavior 
discussed above, ai, not considered explicitly in this paper, could plausibly join n and p eh 
but there is no indication for the a dropping to zero within the range of densities involved. 
How and where they all tend to the mended symmetry limit, if such a limit exists, is not 
clear. This possibility contrasts with the supersymmetric QCD scenario (for p) of [5]. 
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APPENDIX 


A Density-Independence of the Pion Tensor Force 


In Section 14.21 - and in all previous works on tensor forces implemented with IDDs - the 
density dependence of the pion tensor force was ignored, arguing that it is protected by 
chiral symmetry. Here we show explicitly that the argument is correct. 

As shown in Section T4.21 the pion tensor depends on density only via m*. To impose the 
scaling of m* consistently, we take in R-I 
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7T 
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1 + 0.13 * n/no. 


(A.l) 


In R-II, the pion mass must drop fast since the bilinear quark condensate tends to zero, so it 
is reasonable to take it to decease rapidly and vanishing at the VM fixed point. Thus 
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The result is shown in Fig. [T^J We find that the pion tensor is more or less independent of 
density both in R-I and in R-II although m* depends on density, where we used Fermi-Dirac 
distribution function as 
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to make m* be continuous at n = n\ji- Thus taking the pion tensor density-independent in 
doing the Vi ow k calculation in Section [6] is justified. 



Figure 12: (r) = Vf (r) (ti r 2 S 12 ) 1 with n \/2 = 2no- 


31 















B One-Parameter Description of R-I 


In Section [HI we have shown with a small fine-tuning for the constant c/ consistent with 
the expected relation a(N,a) < ci(p,u), nuclear matter can be reproduced well within the 
empirical error bars. Suppose one sticks to the on e-cj strategy and asks how well nuclear 
matter can be reproduced. This has been checked for the range of cj = (0.13 — 0.20). We 
see from Fig. [13] that taking a universal scaling parameter within a narrow range of the cj 
values fails to reproduce Nature. This may be due to different 1/N C corrections contributing 
to the ci coefficients and clearly indicates the extremely fine-tuned nature of the ground state 
property of nuclear matter. 



n ln 0 

Figure 13: Ground-state energy Eo(n,0) of symmetric nuclear matter with one parameter cj 
in the range 0.13 — 0.2. 


C Fate of Hidden Local U (2) Symmetry for (p, co) in Region-II 

In Section 14. 1.21 while local U( 2) symmetry was assumed to hold in Region-I (i.e., n < ni/ 2 ), 
we suggested that it should break down in R-II. There is no known theoretical argument 
either for or against this suggestion. Here we show an unequivocal indication from Nature 
that at least within the present framework the symmetry should indeed break down in the 
density regime n ^ 2n$. The argument is based on the assumption that there is vector 
manifestation (VM) fixed point nyM ~ n c ~ (6 — 7)uq at which the p mass vanishes (in the 
chiral limit). 

Now let us suppose that the U( 2) symmetry holds in R-II. This would imply that near the 
chiral restoration point, the VM would hold for both p and uj. We consider two possibilities: 
One scaling with the same slope in R-II and approaching the same point of VM 

ml/m u « g*/g u « g*/g « (1 -n/n c ) (C.4) 

and the other approaching the VM fixed point with different slopes 

ml/mu « Kg*/g u « k (l -—V (C.5) 

V n c -n l/2 ) 

A drastic simplification is made on both and one should be cautious on the interpretation. 
Nonetheless, the qualitative feature can be taken robust. With all other parameters of Section 
[6]fixed the same, the ground-state energy of symmetric nuclear matter comes out as in Fig. 1141 
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Figure 14: Ground-state energy Eq{u, 0) of symmetric nuclear matter (left panel) and neutron 
matter (right panel) with £7(2) symmetry for (p,u). 


One finds that both the symmetric matter and neutron matter become unstable at n 2 n 0 
for (IC.4|) and n ~ 3no for (|C.5|) . This signals the breakdown. It takes place principally because 
the movement toward the VM fixed point softens the repulsion due to the uj exchange, the 
dropping w-nucleon coupling “winning over” the increase in repulsion coming from the drop¬ 
ping mass. This makes the dilaton-exchange attraction taking over, leading to the collapse. 
Although the over-simplified linear d-scaling must bring about a precocious breakdown, this 
result indicates unequivocally that local U (2) symmetry is untenable at high densities above 
n i/ 2 - We take this as a signal from compact stars that the hidden U(2 ) gauge symmetry 
must be, perhaps badly, broken at high density. 
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